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CHANNEL ESTIMATION FOR OFDM COMMUNICATION SYSTEMS 
INCLUDING IEEE 802.11 A ANB EXTENDED RATE SYSTEMS 

CROSS-REFERENCES TO RELATED APPLICATIONS 
[0001] This application claims the benefit of U.S. Provisional Application No. 60/462,286, 
5 filed April 10, 2003, entitled "Channel Estimation for IEEE 802. 1 1 a and Extended Rates," 
which disclosure is incorporated herein by reference for all purposes. 

[0002] The present disclosure is related to commonly-assigned co-pending U.S. Patent 

Application No. (Attorney Docket No. 02 1245-00141 OUS), filed April 5, 2004, 

entitled "Modified Preamble Structure for IEEE 802.1 la Extensions to Allow for Coexistence 
1 0 and Interoperability Between 802. 1 1 a Devices and Higher Data Rate, MIMO or Otherwise 
Extended Devices," which disclosure is incorporated herein by reference for all purposes. 

BACKGROUND OF THE INVENTION 
[0003] The present invention relates generally to channel estimation in communication 
systems and in particular to channel estimation techniques for wireless communication 
15 systems using multiple subcarriers, including systems compliant with the IEEE 802.1 la 
standard. 

[0004] In a typical wireless communication system in use today, a transmitter and a 
receiver communicate in accordance with a protocol such as the IEEE 802.1 la standard. The 
standard specifies a keying scheme that generally involves grouping bits of data to be 

20 transmitted are grouped and each group of bits is mapped to a symbol (a particular 

combination of phase and amplitude) in a signaling constellation. A number of constellations 
(or keying schemes) are known in the art, including binary phase-shift keying (BPSK), 
quadrature phase-shift keying (QPSK), and quadrature amplitude modulation (QAM) 
constellations. The transmitter includes an encoder that transforms an input data stream into 

25 a sequence of symbols. The transmitter may also insert additional information-carrying 

symbols into the input data stream. For example, protocols such as IEEE 802.1 la specify a 
"preamble" to precede transmitted data. The preamble usually includes parameter values 
related to various features of the signal (e.g., data rate, keying scheme, etc.) as well as a 
known sequence of bits or symbols (referred to herein as a "training sequence") that the 

30 receiver can use for calibration. The symbols are modulated using an appropriate protocol 



(e.g., OFDM (orthogonal frequency division multiplexing) in the case of IEEE 802.1 la) into 
a form suitable for transmission over a channel. The receiver receives the signal, 
demodulates it, and generates an output data stream that ideally is identical to the input data 
stream. 

5 [0005] In practice, the output data stream has a non-zero bit error rate that depends on 
various factors. For example, channel properties such as fading and delay, as well as noise, 
are known to affect the received signal. Receivers in use today generally include circuitry for 
estimating and compensating for channel and noise effects. 

[0006] In typical communication protocols, channel effects are generally estimated by 
10 measuring the receiver's response to a known sequence of input symbols that is inserted into 
the input data stream. For instance, Fig. 1 illustrates a transmission in a system complying 
with the IEEE 802.1 la standard. Each transmission begins with a 20 microsecond (/xs) 
preamble having short training symbols 5, a guard interval LG, long training symbols Z,, and 
a signal field (4 microseconds). The preamble is followed by data. The first eight 
15 microseconds of the preamble comprise ten identical short training symbols S (0.8 (is each) 
that are used for packet detection, automatic gain control and coarse frequency estimation. 
The second eight microseconds comprise two identical long training symbols L (3.2 fis each), 
preceded by a guard interval LG that is the same pattern as the last half (1.6 microseconds) of 
the long training symbol L. The long training symbols can be used for channel estimation, 
20 timing, and fine frequency estimation. The signal field provides information about the data, 
such as the length of the packet, encoding scheme (or data rate), and the like. 

[0007] Fig. 2 shows a long training sequence L\ that is used to generate the signal 
representing the long training symbol in a conventional 802.1 la preamble. This sequence 
represents values (phase shifts) used for each of the 64 subcarriers. As specified in the 

25 standard, the subcarriers span a 20 MHz channel and are spaced apart by 3 12.5 kHz. By 
convention, used here, the first value in the sequence L\ is the value for the DC subcarrier 
(index k= 0), followed by the value for the 1x312.5 kHz subcarrier, then the value for the 
2x312.5 = 625 kHz subcarrier, etc., up to the 32nd value for the 31x312.5 kHz = 9687.5 kHz 
subcarrier. The 33rd value corresponds to the -10 MHz subcarrier, followed by the -(10 

30 MHz - 312.5 kHz) subcarrier, and so on, with the 64th (index k = 63) value corresponding to 
the -312.5 kHz subcarrier. As can be seen from Fig. 2, the DC value and the 28th through 
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38th values (corresponding to subcarriers near the edges of the 20 MHz channel that are not 
used to transmit data) are zero. 

[0008] The output of the transmitter is a long training symbol at a sample rate of 64 
samples/symbol. The samples are obtained by taking a 64-point IFFT (inverse Fast Fourier 
5 transform) of the long training sequence L\. As used herein, a sequence in the frequency 
domain is expressed with uppercase letters (e.g., L[k]), while the corresponding time 
sequence is expressed with lowercase letters (e.g., /[/]). The transformed samples are 
modulated according to the RF channel frequency, converted to analog, and transmitted. 

[0009] The receiver detects the RF signal, converts it to digital (at standard 20 MHz 
10 sampling rate) and demodulates it. Thus, the receiver generates 64 samples during each 
repetition of the long training symbol, for a total of 128 training samples. These training 
samples r[i] (where i = 0, 1, 127) are used to compute channel response according to 
various conventional algorithms. For example, one such algorithm involves averaging the 
corresponding samples from the two repetitions of the long training symbol: 

15 r[i] = (r[i] + r[i + 64])/ 2, i = 0, 1, 63 . (Eq. 1) 

A Fast Fourier Transform (FFT) is performed on the averaged time-domain samples r[i] to 
obtain 64 subcarrier components S[k], k = 0, 1, . . 63, representing the channel amplitude 
and phase for each subcarrier, phase shifted according to the components of the training 
sequence L\. The phase shift due to the training sequence is removed by multiplying each 
20 subcarrier coefficient S[k] by the corresponding training sequence element L\[k], and an IFFT 

is performed on the result to obtain a complex-valued channel impulse response h[i] for each 
subcarrier. 

[0010] Recently, there has been interest in developing multiple input, multiple output 
(MIMO) wireless systems to provide increased data rates. Such systems generally include a 

25 transmitter system having a number (M,) of transmit antennas (transmitters) communicating 
with a receiver system having a number (M r ) of receive antennas (receivers), where M r and 
M t may or may not be equal. In a MIMO communication system, each of the M t transmit 
antennas transmits, at substantially the same time, a symbol representing a different group of 
bits from the input data stream. If each symbol represents B bits, the number of bits 

30 transmitted per symbol period is i?*M,. Thus, MIMO systems are considered desirable as a 
way to increase bandwidth of the communication system. For example, a conventional SISO 
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system compliant with IEEE 802.1 la can transmit up to 54 Mbps (megabits per second). In a 
MIMO system with M t transmit antennas, the data rate can be increased to M/*54 Mbps. 

[0011] MIMO systems introduce additional complexity to the channel estimation that the 
receiver needs to perform. Specifically, each receive antenna receives a signal that is a 
5 combination of signals from the transmit antennas, modified by the channel properties and 
noise. In effect, there is a channel from each transmitter to each receiver, and separate 
channel estimates must be made for each such channel. Thus, the per-carrier channel 
estimate h[i] becomes an M r xM t matrix h[i]. 

[0012] Various techniques for modifying the IEEE 802.1 la preamble for use in MIMO 
10 systems have been suggested. These techniques generally involve changing or repeating the 
preamble to allow separate calibration for each transmit antenna. Some of these techniques, 
such as having each transmit antenna transmit the training sequence while other antennas 
transmit nothing, can lengthen the training period and reduce the data rate. 

[0013] It would be desirable to provide improved channel estimation techniques for use in 
1 5 MIMO-extended 802. 1 1 a and other wireless communication systems. 

BRIEF SUMMARY OF THE INVENTION 
[0014] Embodiments of the present invention provide channel estimation techniques that 
can be used in legacy 802.1 la and extended systems. In some embodiments, channel impulse 
response is adaptively truncated in the time domain during channel estimation, based on 
20 channel power and noise considerations. In other embodiments, channel response tracking 
techniques are implemented to correct for variations in channel response over the 
transmission time of a packet. 

[0015] According to one aspect of the present invention, a method is provided for 
determining channel estimates at a receiver for a wireless communication system using 

25 orthogonal frequency division multiplexing (OFDM) over a plurality of OFDM subcarriers. 
Training signals are received from one or more receive antennas. An estimated channel 
impulse response is computed from the received training signals by reference to a training 
sequence. The estimated channel impulse response is truncated in the time domain to 
improve the signal-to-noise ratio of the channel estimates. In some embodiments, adaptively 

30 truncating the estimated channel impulse response includes computing a cutoff time that 
maximizes a sum of channel power and signal-to-noise ratio; truncating the time domain 
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channel impulse response estimate at the cutoff time; and transforming the truncated time 
domain channel impulse response estimate to a frequency domain channel impulse response 
estimate. 

[0016] According to another aspect of the present invention, a method is provided for 
5 channel estimation for a receiver of a multiple input, multiple output (MIMO) 

communication system wherein signals are transmitted using orthogonal frequency division 
multiplexing (OFDM) over a plurality of OFDM subcarriers. At each of a plurality of receive 
antennas, signals are received from a plurality of transmit antennas, wherein the signal from 
each transmit antenna includes a different subset of the plurality of OFDM subcarriers. The 

1 0 received signals at each receive antenna are transformed to a plurality of impulse coefficients 
for that receive antenna, each impulse coefficient corresponding to a different one of the 
OFDM subcarriers. For each of the receive antennas, a channel impulse response is 
computed for one of the transmit antennas using the impulse coefficients for the subset of the 
OFDM subcarriers transmitted by the one of the transmit antennas. In some embodiments, 

1 5 computing a channel impulse response for one of the transmit antennas includes adaptively 
truncating a time domain channel impulse response estimate to improve the signal-to-noise 
ratio of the channel estimates. 

[0017] According to yet another aspect of the present invention, a method is provided for 
tracking channel variations during receipt of a packet using one or more receive antennas. 
20 An initial channel estimate is determined from training data included in the packet. A 

received symbol in the packet is identified. An input symbol value is estimated using the 
received symbol value and the initial channel estimate. A per-symbol channel estimate is 
derived from the received symbol value and the estimated input symbol value. The initial 
channel estimate is updated using the per-symbol channel estimate. 

25 [0018] The following detailed description together with the accompanying drawings will 
provide a better understanding of the nature and advantages of the present invention. 

BRIEF DESCRIPTION OF THE DRAWINGS 
[0019] Fig. 1 illustrates the structure of a conventional IEEE 802. 1 la packet. 

[0020] Fig. 2 shows a long training sequence LI used to generate a portion of a 
30 conventional IEEE 802. 1 la preamble. 
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[0021] Fig. 3 is a simplified block diagram of a wireless communication system according 
to an embodiment of the present invention, including a transmitter and a receiver. 

[0022] Fig. 4 shows alternative long training sequences suitable for use in an embodiment 
of the present invention. 

5 [0023] Fig. 5 is a flow diagram of a process for obtaining channel estimates for each 

transmitter signal in a MIMO system according to an embodiment of the present invention. 

[0024] Fig. 6 is a flow diagram of a process for obtaining channel estimates with adaptive 
truncation for a SEMO system according to an embodiment of the present invention. 

[0025] Fig. 7 is a flow diagram of an alternative process for obtaining channel estimates 
10 with adaptive truncation for a SIMO system according to an embodiment of the present 
invention. 

[0026] Fig. 8 is a flow diagram of a process for obtaining channel estimates with adaptive 
truncation in a MIMO OFDM system according to an embodiment of the present invention. 

[0027] Fig. 9 is a flow diagram of an alternative process for obtaining channel estimates 
15 with adaptive truncation in a MEMO OFDM system according to an embodiment of the 
present invention. 

[0028] Fig. 10 is a flow diagram of a process for updating channel estimates during receipt 
of a packet according to an embodiment of the present invention. 

DETAILED DESCRIPTION OF THE INVENTION 
20 [0029] Embodiments of the present invention provide improved channel estimation 

techniques that can be implemented in conventional ("legacy") IEEE 802.1 la communication 
systems as well as in various extended rate 802.1 la-based systems and other wireless 
communication systems using OFDM. 

System Overview 

25 [0030] Fig. 3 is a simplified block diagram of a MIMO communication system 300 suitable 
for implementing channel estimation according to an embodiment of the present invention. 
MIMO communication system 300 includes a transmitter system 302 and a receiver system 
304 communicating via a channel 306. Transmitter system 302 includes an encoding block 
312 that receives an input data stream and encodes it into a sequence of symbols (complex 

30 values) to be transmitted, a packet formation block 312 that formats the encoded data into 
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packets for transmission, an IFFT (inverse Fast Fourier Transform) block 316 that applies the 
symbol data as modulation to a number of subcarriers, and a modulator block 318 that 
modulates the IFFT output according to a carrier frequency for transmission via a number 
(M t ) of transmit antennas 320. Herein, multiple instances of like objects are denoted with 
5 reference numbers identifying the object and parenthetical numbers identifying the instance 
where needed. 

[0031] Receiver system 304 includes a demodulator 332 for demodulating a signal received 
via a number (M r ) of receive antennas 330, a channel estimation block 334 that determines 
channel effects from the training data, a decoder 336, and a channel update block 338. 

10 [0032] Communication system 300 is advantageously implemented as an OFDM 

(orthogonal frequency division multiplexing) system, in which signals are split into several 
narrowband channels at different subcarrier frequencies selected to minimize interference 
between the channels. In some embodiments, communication system 300 is operable in 
various modes, one of which is compliant with the IEEE 802.11a standard entitled "Part 11: 

1 5 Wireless LAN Medium Access Control (MAC) and Physical Layer (PHY) specifications, 
High-speed Physical Layer in the 5 GHz Band" incorporated herein by reference. In other 
modes, system 300 may be compliant with various extensions of the IEEE 802.1 la standard, 

such as extensions described in above-referenced Application No. (Attorney Docket 

No. 021245-001410US). 

20 [0033] As shown in Figure 3, input data is provided (generally as a bit stream) to encoder 
312, which encodes the data stream into a sequence of complex numbers (symbols). In some 
embodiments, encoding may include forward error correction (FEC) encoding (which 
involves introducing redundant bits into the data stream), bit interleaving, and mapping 
groups of bits to symbols in a signaling constellation (using, e.g., any of BPSK (binary phase 

25 shift keying), QPSK (quadrature phase-shift keying), 16QAM (quadrature amplitude 

modulation), or 64QAM). Numerous encoding and mapping techniques are known in the art 
and may be implemented in encoder 312. In some embodiments, encoder 312 may also be 
controllable via various parameters to select, e.g., different error correction schemes or 
different signaling constellations. A detailed description of encoding techniques is omitted as 

30 not being critical to understanding the present invention. 

[0034] After mapping, the symbols are provided to packet formation block 314, which 
groups the symbols into packets for transmission. Packet formation block 314 may also add 
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additional symbols into the data stream, such as a conventional IEEE 802.1 la preamble as 
illustrated in Fig. 1 or a modified IEEE 802.1 la preamble, e.g., as described in above- 
referenced Application No. (Attorney Docket No. 02 1245-001 4 10US). 

[0035] Packet formation block 3 1 4 provides groups of symbols to IFFT module 3 1 6, which 
5 uses the group of symbols to modulate a group of OFDM subcarriers. In general, if a group 
of symbols D[k], k = 0, 1, ...,63 is provided, IFFT module 316 generates 64 time-domain 
samples d[i]: 



10 [0036] In an IEEE 802. 1 1 a-compliant embodiment of systems 300, the DC subcarrier and 



1 1 subcarriers near the edges of the 20 MHz channel are not used for data transmission. In 
addition, four of the remaining 52 subcarriers are used to transmit pilot tones. As is known in 
the art, packet formation block 314 may insert appropriate complex values at unused or pilot 
subcarrier locations in each group of symbols D[k] to comply with these aspects of the 



15 standard. In other embodiments, any or all of the subcarriers k may be used to transmit data, 
thereby increasing the data rate. In still other embodiments, even and odd subcarriers (or 
other groupings of subcarriers) may be transformed separately for transmission via different 
transmit antennas 320. 

[0037] Modulator 318 uses the sample values provided by IFFT block 3 16 to modulate an 
20 RF carrier wave and delivers the resulting signal to the M t transmit antennas (transmitters) 
320 for transmission. In some embodiments, each transmit antenna 320 transmits different 
data in parallel, possibly with a cyclic shift in the data from different antennas. In other 
embodiments, each antenna transmits a full set of subcarriers in a different RF channel band. 

[0038] At receiver system 304, the transmitted signal is received, and the reverse of the 
25 transmitter process occurs to reconstruct the input data stream. Thus, demodulator 332 

receives the signals via M r receive antennas 330 (receivers), removes the RF modulation, and 
provides a sequence of samples to the remaining components. In one embodiment, groups of 
samples are provided as an M r -component vector r[z] (with complex components) 
representing the signal sample from each receive antenna 330 at a given time. Each 




(Eq. 2) 



where j = v— 1 . 
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component of r[i] is a linear combination of signals from all of the M t transmitters, with each 
signal being affected by (generally different) channel properties. 

[0039] Channel estimation block 334 is configured to detect a sample stream r[i] 
corresponding to the long training sequence and to use the detected training samples to 
5 generate channel estimates. Recognition of the training sequence may be accomplished using 
a variety of techniques, including mode-sensing techniques described in above-referenced 

Application No. (Attorney Docket No. 02 1245-001 4 10US). The channel estimates 

are advantageously provided as a set of M r xM t matrices M[k] with complex valued 
components that represent the channel effect for a particular transmit-receive pair for a given 

10 subcarrier k. Channel estimation block 334 advantageously implements an improved 
algorithm for generating channel-estimate matrices H[&], such as any of the algorithms 
described below. In some embodiments, channel estimation block 334 may also include 
additional components adapted to estimate other error corrections from the training data, e.g., 
clock and carrier frequency offsets, noise power, and the like. Such aspects of error 

1 5 correction are not critical to the present invention and a detailed description is omitted. 

[0040] Symbol decoder block 336 applies the channel estimates H (as well as any other 
corrections) to the samples r[i] and generates a reconstructed symbol stream x that ideally 
corresponds to the sequence of symbols provided to IFFT block 316. This symbol stream is 
provided to bit stream block 340, which reverses the keying, interleaving, and encoding 
20 processes of encoder block 312 to generate an output bit stream. Symbol decoder block 336 
and bit stream block 340 may be implemented using a variety of techniques, and a detailed 
description is omitted as not being critical to understanding the present invention. 

[0041] In some embodiments, the channel estimates H, the decoded symbols x , and the 
received samples r are provided to a channel update module 338. Channel update module 

25 338 uses the received samples r and the decoded symbols x to detect changes in the channel 
properties that may occur during transmission of a packet. Channel properties may be 
affected by various environmental changes, such as motion of a mobile transmitter system 
302 and/or receiver system 304, operation of other devices in the area, and so on. Channel 
estimation block 334 can be expected to update the channel estimates to reflect such changes 

30 when the next training sequence is received, but in this embodiment, channel estimation 

block 334 does not detect changes that may occur while the data portion of a packet is being 
transmitted. Thus, channel update module 338 advantageously employs adaptive feedback 
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mechanisms (described further below) to compensate for motion or other environmental 
changes. 

[0042] It is to be understood that the MIMO communication system described herein is 
illustrative and that variations and modifications are possible. Any number of transmit 
5 antennas M t and any number of receive antennas M r may be provided, and it is not required 
that M t and M r be equal. In addition, some or all components may be replaced with 
components having a single transmit antenna and/or a single receive antenna. The design of 
various system components, as well as methods used for forward error correction, 
interleaving, and modulation can be varied. The components described herein may be 
10 implemented using hardware, software, or any combination of hardware and/or software 
elements. 

Preamble Generation 

[0043] As described above, packet formation block 314 advantageously prepends a 
preamble to each transmitted packet. In one embodiment, a conventional IEEE 802.1 la 

15 preamble may be used. In other embodiments, a preamble suitable for extended-rate 802.1 la 
communication systems may be used. For example, in a MIMO implementation with two 
transmit antennas (transmitters), each transmitter may transmit a conventional 802.1 la 
preamble while the other transmits nothing. Another option is for each transmitter to transmit 
a subset of the subcarriers (e.g., in a two-transmitter system, one transmitter might transmit 

20 even subcarriers while the other transmits odd subcarriers). 

[0044] Another option, described in above-referenced Application No. (Attorney 

Docket No. 021245-001410US), is to apply a cyclic delay shift on the long training symbol 
and Signal field IFFT outputs prior to applying the guard time extension. For example, 
assume L[k] and D[k] are the 64 subcarrier values for the long training symbol and Signal 

25 field symbol, respectively. For a conventional 802.1 la single transmitter transmission, the 
time samples for the long training symbol are derived by taking the 64-point IFFT of L[k] to 
obtain /[/] and transmitting the samples of /[*]. Thus, with the guard time, the long training 
symbol and guard time are constructed as {/[33:64] /[1:64] /[1:64]}, i.e., the IFFT output is 
repeated twice and the last 32 samples are prepended to form the long training guard interval. 

30 As with the conventional timing, the long training guard interval (32 samples) is twice as 

long as the guard interval for 802. 11a data symbols (16 samples). The signal field is formed 
by {^[49:64] </[l:64]}, where rf[l:64] are the 64 samples of the IFFT of D[k]. 
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[0045] In the case of a two transmitter MEMO device, the first transmitter would transmit 
the long training symbol and signal field as in conventional 802.1 la. The second transmitter 
would apply a cyclic shift such that instead of the IFFT output /[1 :64], it uses the cyclically 
shifted samples fc={/[33:64] /[1:32]} to construct the long training symbol samples as 
5 {/s[33:64] ls[\ :64] ls[l :64]} . For the signal field, the second transmitter uses the shifted 
samples <fc={</[33:64] </[l:32]} to construct the signal field as {^[49:64] *[1:64]}. 

[0046] In a legacy 802. 11a packet, one 3.2 microsecond repetition of the long training 
symbol L is expressed in the time domain as the IFFT of L\[k], where L\[k] contains 64 
subcarrier values, of which 52 are non-zero (as shown in Fig. 2). The time samples are 
10 given as: 



[0047] In an extended mode long training sequence, L\[k] may be modified. For example, 
a modified L[k] can contain more than 52 non-zero subcarriers. Examples of suitable long 
training sequences with additional non-zero subcarriers are shown as sequences L2 and L3 in 

15 Fig. 4, following the same convention for subcarrier order as standard training sequence L\ of 
Fig. 2. L2 is advantageously used in communication systems configured for single input, 
multiple output (SIMO), while Z, 3 is advantageously used in MEMO communication systems. 
Other long training sequences may also be used, including various sequences having one or 
more of the following features: (1) the sequence is formulated such that channel estimation 

20 can be done for multiple transmitters; (2) the sequence has low cross-correlation with the 
standard 802.1 la training sequence; and (3) the sequence is usable in a process for 
determining whether the packet is an 802.1 la packet or an extended-mode packet. Further 
discussion of selection of training sequences may be found in above cross-referenced 



25 [0048] Second, in the case of MEMO transmission, can have a cyclic shift that may be 
different for each transmitter m. The shifted signal / m [f] can be derived from as 
l m [i]=l[(i+64-8 m )%64], where "%" denotes the modulo operator and 8 m is the cyclic delay of 
transmitter m in 20 MHz samples. This expression assumes a 20 MHz sampling rate, such 
that there are 64 samples in a 3.2 microsecond interval. An alternative method of generating 

30 the cyclic shift is to apply a phase ramp rotation to all subcarrier values of L(k) prior to 
calculating the IFFT, e.g.: 




(Eq. 3) 



Application No. 



(Attorney Docket No. 021245-001410US). 
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[0049] A MIMO transmitter system can have two or more transmit antennas (or antenna 
arrays, as the case may be). For a MIMO system with two transmit antennas and two 
different transmit data streams, preferred values for the cyclic delay values 8 m are 0 and 32 
5 samples, respectively. This corresponds to a cyclic delay of 1 .6 microseconds between the 
two transmitters. For three transmitters, 5 m can be 0, 22, and 43 samples, respectively. For 
four transmitters, 8 m can be 0, 16, 32, and 48 samples, respectively. 

Overview of Channel Estimation in a MIMO System 

[0050] At the receiver side, the channel estimates for each transmitter signal can be 
10 estimated by a process such as process 500 shown in Fig. 5. Process 500 begins with 

receiving and (digitally) sampling signals for the long training symbol (step 501). Then, a 
64-point FFT (Fast Fourier Transform) of the received long training symbol samples is taken 
(step 502), as is done for conventional 802.1 la preamble reception. Next, each subcarrier is 
multiplied by known phase shift (Z) values (step 503), and an IFFT of the result is taken to 
15 get a 64-point impulse response estimate (step 504). 

[0051] In one embodiment of a MIMO transmission, these 64 samples contain the cyclic 
shifted impulse responses for all M t transmitters. With that, the receiver can isolate the 
impulse responses for each MIMO transmitter (step 505). For MIMO with two transmit 
streams, this can be done by separating the first 32 samples and last 32 samples. For four 
20 transmit streams, groups of 16 samples can be extracted. 

[0052] From the extracted impulse responses per transmitter, channel estimates can be 
derived (step 506) for all subcarriers by taking a 64-point FFT of each impulse response, 
where zero values are appended to the sample values to provide 64 input values. Where the 
receiver system has multiple receive antennas (or antenna arrays), process 500 may be 
25 performed for each receive antenna, resulting in a complete matrix of channel estimates. 

Channel Estimation in a SEMO System 

[0053] Some extended rate systems are single input, multiple output (SIMO) systems, in 
which one transmit antenna sends data to a number of receive antennas. Channel estimates 
are generated for each receive antenna. In some embodiments, channel impulse response 
30 time is generally shorter than the symbol period. For example, in IEEE 802. 11a systems, the 
symbol period is 4 /is, while the channel impulse response time might be much shorter. 
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Channel estimation block 334 can make use of this fact to adaptively truncate the impulse 
response, thereby improving the signal-to-noise ratio of the channel estimates. 

[0054] Fig. 6 is a flow diagram of a process 600 for generating channel estimates using 
adaptive truncation in a SIMO system in which one transmit antenna sends a data signal that 
5 is received by a number M r of receive antennas. 

[0055] At step 602, samples for the long training symbol are received. In an 802.1 la- 
compliant system, samples are received at 20 MHz, and each sample is provided as a 
complex number representing amplitude and phase. Since the symbol is transmitted twice, 
with a duration of 3.2 /xs each time, 128 samples r n [i] are received for each receive antenna n 
10 (n = 0, ...,M r -l). 

[0056] At step 604, corresponding samples from the first and second transmission of the 
long training symbol are averaged: 

s.[i] - SK±ll±aE±I±gl , , = 0, .., 63. (Eq. 5) 

2 

Here, ris an estimated starting sample for the channel estimation, which may be, e.g., about 
15 16 to 32 samples after the start of the long training symbol. 

[0057] At step 606, a 64-point FFT of s[i] is computed to provide 64 subcarrier coefficients 
S n [k],k=0, 1, ...,63: 

S„[k] = f>„ Wexpf-;^ . (Eq. 6) 



i=0 



[0058] At step 608, the subcarrier coefficients S n [k] are multiplied by the respective long 
20 training sequence values L[k] (for the long training sequence in use for a given embodiment, 
e.g., conventional sequence L\ of Fig. 2 or extended sequence L 2 of Fig.4 ) to determine 

impulse coefficients S° n [k] : 

S° n [k] = L[k]S n [k], £ = 0,...,63. (Eq. 7) 

[0059] Because L[0] is zero, the impulse coefficient at DC is "missing." At step 610, an 
25 impulse coefficient at DC is interpolated from the other impulse coefficients: 



32 



(Eq. 8) 
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[0060] At step 612, a 64-point IFFT is computed for the spectral coefficients S°[k] , 
resulting in a channel impulse response estimate: 

/U0 = IX[*]exp y25* . (Eq.9) 

[0061] To increase the signal-to-noise ratio (SNR) of the channel estimates, the impulse 
5 response of Eq. 9 is advantageously cut off (truncated) in time at a point that optimizes a 
tradeoff between channel power and noise. More specifically, in this embodiment, at step 
613, a channel power function P h [t] is defined for a time index f = 0, ...» 63 as: 

= S 5* !<,[(/ + 0 mod 64] 2 , (Eq. 10) 

where mod is the modulo operator. This function approximates the channel power received 
10 during a relatively short time interval beginning at time t (measured in samples). In other 

embodiments, the sum over 8 samples may be replaced with a sum over a different number of 
samples, and the modulo operator may be replaced with a cutoff at i + t = 64. 

[0062] At step 614, a channel noise power estimate cr 2 that represents the amount of noise 
included in the channel power loss function P h [t] is generated. A variety of techniques can 

15 be used to generate this estimate; in one embodiment, 

<r 2 =minP A [f]. (Eq. 11) 

r 

Eq. 1 1 is a reliable noise estimate to the extent that P h [t] reflects a channel impulse response 

superimposed on a substantially constant background of noise, and to the extent that at some 
time index t, the channel response is dominated by the noise. It will be appreciated that other 
20 estimates of channel noise may also be used. 

[0063] At step 615, a cutoff time t c is computed as: 

^=argmin(PJ^] + ^ 2 ),/ = ^ 0 ,...,63, (Eq. 12) 

k 

where P h [t] is given by Eq. 10, cr 2 is the channel noise power estimate from step 614, and 

,! arg min" indicates that t c is the value of t for which the sum is minimized. A minimum 
25 cutoff time to is advantageously predefined based on a minimum impulse response length 
determined by all filters at both transmit and receive ends of the channel. In one 
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embodiment, t 0 = 12. More generally, a variety of methods for determining a minimum 
impulse response length are known in the art and may be used to define to. 

[0064] It should be noted that Eq. 12 represents a tradeoff between keeping more or fewer 
of the time samples for the channel impulse response. As fewer samples are kept, the noise is 
5 reduced but more of the useful channel response information is lost. Other algorithms for 
evaluating this tradeoff to select a cutoff time may be substituted for that described herein. 

[0065] At step 616, the 64-point FFT of the truncated impulse responses is computed as: 



63 



# B °[*] = 2> n C Mexp 



J 



2nik 
64 



(Eq. 13) 



where: 



10 h c n \i} = \ (Eq. 14) 

0, i>t c 

is the truncated channel impulse response estimate. 

[0066] The FFT outputs of Eq. 13 are channel estimates H 0 [k] that may be used for 

channel correction of received data signals. Various channel correction techniques are 
known in the art and may be employed in conjunction with channel estimation process 600. 

15 [0067] Fig. 7 is a flow diagram of an alternative process 900 for channel estimation with 
adaptive truncation in a SIMO OFDM system. This process is generally similar to process 
600 described above except that the DC offset estimation is done somewhat differently and 
the computation of the cutoff time is modified. In addition, process 900 introduces a timing 
refinement. 

20 [0068] More specifically, steps 902, 904, 906, and 908 correspond to steps 602, 604, 606, 
and 608 of process 600 described above and result in impulse coefficients S*[k] in 

accordance with Eqs. 5-7 given above. At step 910, a 64-point EFFT of the impulse 
coefficients is computed to obtain channel impulse response estimates as: 

*=o V 64 ) 
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It should be noted that Eq. 15 differs from Eq. 9 above in that the DC coefficient S%[k] has 

not been interpolated to a nonzero value. Instead, in process 900, the DC offset is corrected 
in the time domain. Specifically, at step 912, the DC offset for each receiver is estimated by 
averaging the tail of the channel impulse response of that receiver over 16 samples: 



5 ^-ifliP]. < E <1- 16 ) 

lO i=40 

where the sample window (40 to 55 in Eq. 16) may be optimized for a particular embodiment 
as desired. The sample window is advantageously chosen to be near the end of the sampling 
interval but may be shifted away from the end somewhat, e.g., to avoid effects of any timing 
inaccuracy in choosing the starting sample. In addition, the size of the sample window may 
10 be varied as desired. 

[0069] At step 914, the channel impulse response estimates h n [i] are corrected for the DC 
offset: 



Or « 



h° n [i] = K[i]-x n . (Eq. 17) 

[0070] Next, the timing estimate rthat was used to select the sample window (see Eq. 5) is 
15 refined based on the channel impulse response. At step 916, a channel power function is 
computed as: 

7 M r -\ ^ 2 

^] = ZZ A. 0 [(i + /)mod64] , (Eq. 18) 

i=0 n=0 

similarly to step 613 of process 600. At step 918, the timing correction r' is estimated as: 

r' = argmaxP A |>]> * = 0,...,63; (Eq. 19) 

t 

20 that is, r' is the value of t for which P h [t] is maximized. This timing correction is used to 
adjust the initial sample time from r to r'. Accordingly, at step 920, the impulse responses 
and the power estimate are cyclically shifted by — r' samples: 



ftp] = ft K< + O mod 64] , (Eq. 20) 



and 



25 P A r [/] = PJ(/ + r')mod64]. (Eq. 21) 
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[0071] Next, at step 922, a cutoff time t c for the impulse response is determined. The 
cutoff time can be determined according to Eq. 12 above. In an alternative embodiment, the 
noise power estimate is obtained as: 



cr n 2 =min(/> r [/])/4, (Eq. 22) 



and the cutoff time t c is given by: 



r c =min{P;(/)<2a n 2 }, / = 8,...,31, (Eq. 23) 

i 

that is, the smallest time index i for which P*[i] is less than twice the noise power estimate. 

The minimum cutoff time (8 in Eq. 23) is selected to prevent unrealistically small cutoff 
values from being used in channels with large noise or interference levels. As in process 600, 
this minimum time may be selected based on the known filtering properties of the transmit 
and/or receive antennas. 

[0072] At step 924, the time-corrected channel impulse responses h T n (i) are truncated. In 

some embodiments (e.g., where the transmitter is a legacy transmitter or other type of 
transmitter that does not use all subcarriers), abrupt truncation in the time domain can give 
rise to Gibbs* phenomena in the frequency domain. To mitigate this effect, the truncation 
may be done using a windowing function w(i) defined as: 



w(/) = 



1 i<t c -4 

- (1 + cos(/r(z - t c - 4) / 8)) t c - 4 < i < t c + 4 
2 

0 f„+4</<58 

-(l + cos(^(64-/)/8)) 58 <i 

k. 2 



(Eq. 24) 



The channel impulse response is then given by: 



It will be appreciated that other windowing functions may be substituted. For instance, in 
embodiments where all subcarriers are used, the windowing function w(i) = 1 for all i may be 
used. 

[0073] At step 926, final channel estimates (in the frequency domain) are computed as a 
64-point FFT of the truncated channel impulse responses of Eq. 25: 



17 



. 2mk ) 

f 

64 J 



(Eq. 26) 



[0074] In some embodiments, at step 928, the channel estimates for subcarriers near the 
channel edges may be replaced by untruncated channel estimates. For example: 



5 may be used. 

[0075] It will be appreciated that the channel estimation processes described herein are 
illustrative and that variations and modifications are possible. Steps described as sequential 
may be executed in parallel, order of steps may be varied, and steps may be modified or 
combined. Process 600 or 900 may also be implemented in a conventional (legacy) 802.1 la 

10 system including single receiver systems. In addition, various steps or equations of process 
600 may be combined with or substituted for corresponding steps or equations of process 
900, or vice versa. For example, the windowing function of Eq. 24 may be used to truncate 
channel estimates in a process generally conforming to process 600. These processes may 
also be extended to MIMO systems, e.g., where different transmit antennas transmit the 

15 training sequence at different times. 

[0076] Any desired training sequence L[k] may be used, including L 2 or L3 shown in Fig. 4. 
For example, process 600 or 900 may be implemented in a legacy 802.1 la communication 
system, e.g., with a single transmitter and a single receiver using the standard training 
sequence L\[k] of Fig. 2; in this instance, the window function of Eq. 24 is advantageously 

20 applied. As another example, in an extended rate system where all subcarriers are used, 

training sequence L 2 [k] of Fig. 4 might be used, and the window function of Eq. 24 might not 
be applied. In some embodiments, where all subcarriers except at DC are used to transmit 
data, a modified training sequence such as L 2 [k] or L^[k] (Fig. 4) that includes nonzero 
elements for all subcarriers may be used. It should also be noted that a modified training 

25 sequence such as L 2 [k] or L 3 [k] may indicate to the receiver that the transmitter is using an 
extended rate transmission protocol, as described in above-referenced Application No. 



H° n [k] = S° n [k]cxp -j 



. 2nkt c \ 
f - 

64 ) 



£ = 26,27,39,40 



(Eq. 27) 



(Attorney Docket No. 021245-001410US). 
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Channel Estimates for MIMO OFDM Systems 

[0077] For MIMO systems with multiple independent transmitters, a channel estimate is 
advantageously computed for each transmitter-receiver pair. For example, in an embodiment 
of system 300 (Fig. 3) with two transmit antennas 320 and two receive antennas 330, four 
5 channels are present: (1) from antenna 320(1) to antenna 330(1); (2) from antenna 320(2) to 
antenna 330(1); (3) from antenna 320(2) to antenna 330(1); and (4) from antenna 320(2) to 
antenna 330(2). Process 600 described above may be adapted to generate a channel estimate 
for each of these channels. 

[0078] Fig. 8 is a flow diagram of a process 700 for computing channel estimates for each 
10 of a number of receive antennas in a MIMO system. In this example, there are two transmit 
antennas (m = 0 and m = 1), and each transmit antenna transmits only half the subcarriers for 
the training sequence. For example, transmit antenna m — 0 may transmit the even 
subcarriers while transmit antenna m = 1 transmits the odd subcarriers. Process 700 involves 
splitting the subcarrier coefficients S n [k] for each receive antenna n into two subsets S nm [k] 
1 5 and deriving a separate set of channel estimates for all subcarriers from each subset. 

[0079] More specifically, at step 702 (which may be generally similar to step 602 of 
process 600 described above), samples for the long training symbol are received for each 
receive antenna. At step 704, for each receive antenna «, samples from the first and second 
half of the long training period are averaged, similarly to step 604: 

20 •?„!>] = — — - -,i = 0, ...,63. (Eq. 28) 

2 

[0080] At step 706, a 64-point FFT of s n is computed for each receive antenna: 

S n [k] = £ s n [i] exp [-J . (Eq. 29) 

[0081] At step 708, the 64 subcarrier coefficients S„[k] for each receive antenna n are split 
into two subsets S nm (where m = 0, 1 denotes the transmit antenna), each including 32 
25 coefficients In one embodiment, the even subcarriers (k = 2p for integer p) are associated 
with transmit antenna 0 and odd subcarriers {k = 2/H-l for integer p) are associated with 
transmit antenna 1 . 

[0082] At step 710, the subcarrier coefficients are multiplied with the known long training 
sequence values L[k] to determine impulse coefficients S° m [p] : 
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S° n0 [p] = L[2p]S tt [2p] 

(Eq- 30) 

S° al [p] = L[2p + l]S a [2p + l] 



for p = 0, 1 , . . . , 3 1 . In embodiments where L[0] is zero, S° 0 [0] becomes a "missing" impulse 
coefficient. At step 712, the missing impulse coefficient is supplied by interpolation, using: 

Si[0] = IS>]expf-;^l. (Eq. 31) 

Thus, 32 impulse coefficients are available for each transmit-receive pair. 

[0083] At step 714, a 32-point IFFT of each set of impulse coefficients S° m [p] is computed 
according to: 

This provides a channel impulse response estimate for each receiver-transmitter pair nm. 

[0084] To increase the SNR of the channel estimates, the impulse responses are 
advantageously cut off (truncated) in time at a point that optimizes a tradeoff between 
channel power and noise. Specifically, at step 715, a channel power loss function Pt\t\ is 
defined for a time index t = 0, . . ., 63 as: 

7 M r -lM ( -\ 2 

JU'] = Z£Z£K' + 0"KMi32] , (Eq.33) 

/=0 n=0 m=0 

where M r is the number of receive antennas and M t is the number of transmit antennas. This 
function approximates the channel power received (from all transmitters) during a relatively 
short time interval beginning at time t (measured in samples). In other embodiments, the sum 
over 8 samples may be replaced with a sum over a different number of samples, and the 
modulo operator may be replaced with a cutoff at i + t = 32. 

[0085] At step 716, similarly to step 614, a channel noise power estimate a 2 that represents 
the amount of noise included in the channel power function P h [t] is generated. A variety of 
techniques can be used to generate this estimate; in one embodiment, 

a 2 = nun />[*]. (Eq. 34) 
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r 



Like Eq. 1 1 above, Eq. 34 is a reliable noise estimate to the extent that Ph[t] reflects a channel 
impulse response superimposed on a substantially constant background of noise, and to the 
extent that at some time index /, the channel response is dominated by the noise. Other 
estimates of channel noise may also be used. 

5 [0086] At step 717, a cutoff time t c is computed as: 

t c = min(P h [t] + t<T 2 ) 9 t = t o9 ... 9 63 (Eq. 35) 

where Ph[t] is given by Eq. 33 and a 2 is the channel noise power estimate from step 716. As 
in process 600, a minimum cutoff time to is advantageously predefined based on a minimum 
impulse response length determined by all filters at both transmit and receive ends of the 
10 channel. In one embodiment, to = 12. More generally, a variety of methods for determining a 
minimum impulse response length are known in the art and may be used to define t 0 . 

[0087] At step 718, zeros are appended to the truncated channel impulse response to 
provide a total of 64 samples, and the (64-point) FFT is computed as: 



7=o V 64 j 



1 5 where 



h\i\ i < t r 

£[«> r C - (Eq-37) 

0, L < i < 63 

The FFT outputs of Eq. 36 are the channel estimates H° nm [k] for transmitter m, receiver n, and 

subcarrier k\ these channel estimates may be used for channel correction of received data 
signals. Various channel correction techniques for MIMO systems are known in the art and 
20 may be employed in conjunction with channel estimation process 700. 

[0088] It should be noted that process 700 results in channel estimates for all 64 subcarriers 
k for each receive-transmit pair nm. Thus, during transmission of the data portion of a 
packet, each transmit antennas may use all of the subcarriers, or any subset thereof, including 
overlapping subsets, and appropriate channel corrections may be applied. 

25 [0089] Fig. 9 is a flow diagram of an alternative process 1000 for channel estimation in a 
MIMO extended rate system. This process is generally similar to process 700 described 
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above except that the DC offset estimation is done somewhat differently and the computation 
of the cutoff time is modified. In addition, process 1000 introduces a timing refinement. 

[0090] More specifically, steps 1002, 1004, 1006, 1008 and 1010 correspond to steps 702, 
704, 706, 708 and 710 of process 700 described above and result in impulse coefficients 

5 S° m [k] in accordance with Eqs. 15-17 given above. At step 1012, a 32-point BFFT of the 

impulse coefficients for each transmit-receive pair is computed to obtain channel impulse 
response estimates as: 



31 f 2m(2k + m)^ 



KM - 5XM«P J V/ ' ' CBq. 38) 

*=o V o4 J 

Alternatively, a 64-point EFFT may be used; in that case, S° m [k] is set to zero for odd 
10 subcarriers for transmit antenna m = 0 and for even subcarriers for transmit antenna m = 1. 

[0091] It should be noted that Eq. 38 differs from Eq. 32 above in that the DC coefficient 
S° 0 [k] has not been interpolated to a nonzero value. Instead, in process 1000, the DC offset 

is corrected in the time domain. Specifically, at step 1014, the DC offset for each receiver is 
estimated by averaging the tail of the channel impulse response of that receiver over 4 
15 samples: 

1 28 

*-=7Z$.oM (Eq-39) 

where the sample window (25 to 28 in Eq. 39) may be optimized for a particular embodiment 
as desired. The sample window is advantageously chosen to be near the end of the sampling 
interval but may be shifted away from the end somewhat, e.g., to avoid effects of any timing 
20 inaccuracy in choosing the starting sample. In addition, the size of the sample window may 
be varied as desired. 

[0092] At step 1016, the channel impulse response estimates h n0 [i] are corrected for the 
DC offset: 

&M = 4oM-*.o ( Ec i- 40 ) 

25 It should be noted that DC offset correction is applied only for a transmit antenna whose 
spectral coefficients include a DC term, in this case transmit antenna m = 0. Thus, in this 

example, h al [i] = h n] [i] . 
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[0093] Next, the timing estimate rthat was used to select the sample window (see Eq. 28) 
is refined based on the channel impulse response. At step 1018, a channel power function is 
computed as: 

m = i Z 2|^l[0' + <)mod32]| 2 , (Eq. 41) 

1=0 «=0 m=0 

5 similarly to Eq. 33 above. At step 1020, the timing correction r' is estimated as: 

r' = arg max />[*], t = 0,...,32 , (Eq. 42) 

similarly to Eq. 19 above. This timing correction is used to adjust the initial sample time 
from r to r'. Accordingly, at step 1022, the impulse responses and the power estimate are 
cyclically shifted by -r' samples: 

1 0 VL [i] = hi [0* + r ') mod 32] , (Eq. 43) 



and 



K [/] = P h W + *') mod 32] . (Eq. 44) 



[0094] Next, at step 1024, a cutoff time t c is determined. The cutoff time can be 
determined according to Eq. 34 above. In an alternative embodiment, the noise power 
1 5 estimate is obtained as: 



C7 2 =min(P A r [/])/4, (Eq.45) 



and the cutoff time t c is given by: 



t c = min{P;(/) < 2a- 2 }, i = 8, 31 , (Eq. 46) 

that is, the smallest time index i for which P£[i] is less than twice the noise power estimate. 

20 The minimum cutoff time (8 in Eq. 46) is selected to prevent unrealistically small cutoff 

values from being used in channels with large noise or interference levels. As in process 700, 
this minimum time may be selected based on the known filtering properties of the transmit 
and/or receive antennas. 

[0095] At step 1026, the final channel estimates (in the frequency domain) are computed as 
25 a 64-point FFT of the truncated channel impulse responses, after appending 32 zeroes: 
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= exp -j 



63 / 



. 2nik 

i 

64 



(Eq. 47) 



where 



r 



0, t c < i < 63 




(Eq. 48) 



[0096] In this embodiment, all subcarriers are used to transmit data (including training data 
5 of a sequence such as Li or L 2 of Fig. 4), and a windowing function is not needed for optimal 
results. In other embodiments, a windowing function (e.g., similar to Eq. 24 above) may be 
introduced if desired. 

[0097] In some embodiments, the signal field of the packet (see, e.g., Fig. 1) may be 
transmitted in SIMO mode. In that case, the MIMO channel matrix H Q nm [k] may be used to 
10 obtain a (vector) set of SIMO channel coefficients for proper demodulation of that portion of 
the packet. Specifically, at step 1028, a SIMO channel vector H° n [k] is computed as: 



where 9 = (2^-/64) , k is the OFDM subcarrier index, and 5 is the delay (in samples) applied 
to the OFDM symbol transmitted by the second (m = 1) transmitter to achieve transmitter 
1 5 diversity. 

[0098] It will be appreciated that processes 700 and 1000 are illustrative and that variations 
and modifications are possible. Steps described as sequential may be executed in parallel, 
order of steps may be varied, and steps may be modified or combined. The computations 
may be performed sequentially or in parallel for each receive antenna n. Any desired training 
20 sequence may be used, including sequence L 2 or L 3 shown in Fig. 4. 

[0099] Process 700 or 1000 may also be adapted to a larger number of transmit antennas m 
by dividing the 64 available subcarriers (or other number of subcarriers) as equally as 
possible among the transmit antennas during the training period. It will be appreciated that 
for a sufficiently large number M t of transmit antennas, 64/M/ subcarriers might not provide 
25 enough information for a reliable channel estimation. In such embodiments, the training 

period may be extended by repeating the long training symbol (e.g., using different subsets of 
the transmit antennas for each repetition), or different antennas may transmit uncorrelated 



HSk] = Hl[k] + ei 8ke Hlm, 



(Eq. 49) 



training sequences using the same subcarriers, or other techniques may be employed to 
improve reliability of the channel estimates. 

Channel Variation Tracking 

[0100] In some embodiments, channel tracking may be implemented to improve the 
5 robustness to channel variations during a packet caused by motion of the user or other 

changes in the environment as mentioned above. For example, as mentioned above, channel 
update module 338 of receiver system 304 (Fig. 3) may provide channel tracking and 
adaptive updating during reception of a packet. 

[0101] Fig. 10 is a flow diagram of a process 800 for updating channel estimates during 
10 receipt of a packet. The process is described in the context of an 802.1 la packet and a single- 
antenna receiver; extensions to other systems will be apparent. 

[0102] At step 802, during the training period, an initial channel estimate H Q [k] is 

computed for each subcarrier k. Step 802 may include executing process 600 or 700 
described above, or another channel estimation process. 

15 [0103] At step 804, during receipt of the data portion of the packet, a sequence of samples 
r[i] corresponding to a symbol period is received. At step 806, this sequence of samples r[i] 

is used to estimate a transmitted symbol X[k] for each subcarrier k. Step 806 may be 

implemented using conventional decoding technologies and advantageously includes 

correcting for channel effects using initial channel estimate H 0 [k] . 

20 [0104] At step 808, the samples r[i] for the symbol period are transformed to a set of 

received subcarrier coefficients R[k]; these coefficients will generally reflect the transmitted 

symbols as modified by channel effects. At step 810, a "current" channel estimate H e [k] is 

computed as: 

H e [k] = R[k]X~ l [k]. (Eq.50) 

25 [0105] At step 812, the current channel estimate can be used to update the initial channel 
estimate H 0 [k] , thereby providing an updated channel estimate H'[k] . Various filters may 
be used to generate the updated channel estimate. For example, a first order filter: 

H'[k] = (1 - a)H 0 [k] + aH e [k], 0 < a < 1 (Eq. 51) 
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may be applied. Since H 0 [k] and H e [k] will be similar to the extent that the channel 
properties have not changed since the training period, Eq. 5 1 may change the channel 
estimate only slightly (or not at all). The averaging constant a may be varied as desired; in 
one embodiment a value of a = 0.125 is advantageously used. More generally, the value of a 
5 may be optimized for a particular embodiment based on tradeoffs between sensitivity to 
actual changes in the channel properties and undesired responsiveness to transient 
fluctuations. 

[0106] In other embodiments, updated channel estimate H f [k] may be computed using 
other finite or infinite impulse response filters, e.g., filters that average H e [k] over a number 
10 of symbols to further reduce sensitivity to transient fluctuations in channel response. 

[0107] At step 814, the updated channel estimate H'[k] is used to decode subsequent 
packets. Steps 804 to 812 may be performed again using updated channel estimate H'[k] in 
place of initial channel estimate H 0 [k] , thereby allowing the receiver to respond to further 
changes in channel properties. 

1 5 [0108] For systems with multiple receive antennas, process 800 may be performed 
(sequentially or in parallel) for each receive antenna. 

[0109] For extended rate systems implemented using MEMO, where different data symbols 
are transmitted simultaneously from different transmit antennas, the above technique can be 
used for every group of M t symbols, where M t is the number of transmit antennas. In one 

20 such embodiment, M 0 [k] is an M r xM t matrix with the channel coefficients for subcarrier k, 

X[k] is an M t xM r matrix where each of the M r columns contains the M t estimated symbols 
(one per transmit antenna) transmitted for subcarrier k, and K[k] is a vector with M r complex 
values for all receivers for subcarrier k. Corresponding matrix operations can be substituted 
for the scalar operations in Eqs. 50 and 51. 

25 [0110] While the invention has been described with respect to specific embodiments, one 
skilled in the art will recognize that numerous modifications are possible. For instance, the 
techniques described herein may be adapted to various OFDM communication systems and 
protocols and are not limited to any particular number or spacing of subcarriers. For 
example, FFTs and IFFTs with different numbers of points could be used, or time/frequency 
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transform algorithms other than FFT, such as discrete cosine transforms or the like, could be 
substituted. The techniques may be implemented for use in systems with any number of 
transmit antennas and any number of receive antennas, some or all of which might be mobile. 
Devices implementing the channel estimation techniques described herein may also be 
5 configured for use in legacy 802.1 la networks as well as extended rate networks and may be 
components of a "mixed" network in which legacy and extended rate devices coexist or 
interoperate. 

[0111] The various functional blocks described herein may be implemented using hardware 
(e.g., suitable integrated circuit components), software (e.g., program instructions) executing 
10 on appropriate processing devices, or any combination thereof. In some embodiments, a 
single device may be implemented as a transceiver including both transmitter system and 
receiver system elements. In such devices, hardware and/or software components may be 
configured to perform both transmit and receive functions. 

[0112] Thus, although the invention has been described with respect to specific 
1 5 embodiments, it will be appreciated that the invention is intended to cover all modifications 
and equivalents within the scope of the following claims. 
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